arXiv: 1501.01928v2 [quant-ph] 18 Aug 2015 


Geometric Phase Atom Optics and Interferometry 

B. ZygelmarJ^ 

Department of Physics and Astronomy, University of Nevada, Las Vegas, Las Vegas NV 89154 

(Dated: August 19, 2015) 

We illustrate how geometric gauge forces and topological phase effects emerge in atomic and 
molecular systems without employing assumptions that rely on adiabaticty. We show how geo¬ 
metric magnetism may be harnessed to engineer novel quantum devices including a velocity sieve, 
a component in mass spectrometers, for neutral atoms. We introduce and outline a possible ex¬ 
perimental setup that demonstrates topological interferometry for neutral spin 1/2 systems. For 
that 2-level system, we study the transition from Abelian to non-Abelian behavior and explore its 
relation to the molecular Aharonov-Bohm (MAB) effect. 


I. INTRODUCTION 

Berry’s phase[T], a realization of a non-integrable phase 
factor [5J13], plays an important role in describing adia¬ 
batic quantum evolution in a semi-classical setting. It 
is also known[4H8| that effective gauge potentials, that 
give rise to such phases, emerge in fully quantal systems 
in which true adiabaticity is ill-defined. In addition to 
generating phase holonomies[3], they may lead to effec¬ 
tive Lorentz-like forces [b] H] |9| acting on the quantum 
system. 

Today, applications of geometric gauge forces in cold 
many-body systems [131 El is an active and topical 
area of research. Dressing atoms using lasers [T2|, 
researchers [131 114| have been able to engineer Lorentz- 
like forces, an effect sometimes called geometric [31 ITS] . 
synthetic El or artificial magnetism [10], in ensembles of 
cold atoms. It is hoped that the latter may allow re¬ 
alization of novel quantum Hall physics in a quantum 
degenerate gas. 

Another possible application of geometric magnetism 
is in the manipulation of individual neutral atoms or 
neutrons [T5HT7| . In Refs. El E] we proposed how 
gauge forces could be exploited to construct “magnetic 
lenses” for neutral matter. In this paper we elaborate on 
those observations by introducing two novel and striking 
illustrations of the latter. In the first example we illus¬ 
trate how this phenomenon allows realization of a veloc¬ 
ity sieve, a component in mass spectrometers, for beams 
of neutral particles. In another example we consider a 
quantum mechanical analog El [T31 E] of a field the¬ 
oretical model El to demonstrate topological quantum 
interferometry for neutral spin-1/2 systems. It’s labo¬ 
ratory realization could have applications in topological 
quantum computing protocols. 

With some exceptions e.g. E1II2112Q], most theoreti¬ 
cal studies, by necessity, have relied on some form of the 
adiabatic assumption, i.e. the Born-Oppenheimer (BO) 
approximation. Gauge potentials are explicit in the adia¬ 
batic picture but it also known El [22] that those effective 
(non-Abelian) gauge potentials describe a pure gauge[3]. 
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Consider the following Schroedinger equations 

- = ( 1 ) 

and 

- 2 ^ ( Vh - *A) V' + V{R)^' = ih^ (2) 

for amplitudes if>, ip' respectively. We take them to be 
n dimensional column vectors and V {R) is a local n x n 
matrix potential. The gauge potential 

A = iU\R)VRU{R) (3) 

where U{R) is a differentiable, single-valued, unitary 
n X n matrix. According to definition (§ A is a pure 
gauge mm We pose the question; do the Schroedinger 
equations (Q and ([^ describe the same physics? If 
n = I the answer is in the affirmative since then U (R) = 
exp(—zH) and 

A = i exp(zH) Vr exp(—zH) = VijD (4) 

where H is a single-valued function of R. Therefore, 
■0' = exp(zH) Ip and the physical content of Eqs ([^ and 
(U are identical [33]) as the eigenvalues of physical op¬ 
erators, which transform in a covarient manner, are in- 
varient under this gauge transformation. The same is 
true for the non-Abelian case where rz > 1, provided that 
[U,V] =0 since then p' = U{R)p. So for these cases we 
find that the minimal coupling of a pure gauge potential 
is fully equivalent to a description where A = 0. However 
this is not the case if [A, V] 7^ 0, as the amplitudes p, 
and p' are no longer related by the (single-valued) gauge 
transformation U. The conclusion follows from the fact 
that [U{R), V] 7^ 0 if the above inequality holds (see also 
Appendix C). The implications of this observation will 
become apparent in the discussion below. 

Are non-trivial gauge forces an artifact of the adia¬ 
batic approximation[20] ? Are singularities in the adi¬ 
abatic Hamiltonian solely responsible for the emergence 
of topological phases? Several interpretations [33H1S] for 
the origins of geometric gauge forces have been advanced. 
However, compelling examples that offer fully quantum 
solutions to systems in which such forces arise in the 
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adiabatic limit have largely been unavailable and, there¬ 
fore, predictions are limited by the validity of assump¬ 
tions based on adiabaticity. 

In order to address questions and dehciencies in theo¬ 
retical approaches that assume adiabaticity, we nulls] 
introduced a wave packet propagation scheme that does 
not require the assumptions and restrictions imposed by 
it. The resulting time dependent solutions for the sys¬ 
tems are exact, within the bounds imposed by numerical 
error, and they are not compromised by issues relating 
to the robustness of the adiabatic approximation. This 
allows us to make definitive verdicts on the fidelity of 
predictions informed by the latter, and the gauge theory 
interpretation that follows from it. 

In the discussion below we hrst present the general the¬ 
oretical framework for the specific cases considered. We 
introduce a 2-level system whose adiabatic Hamiltonian 
is defined in such a way that there is a cancellation of 
the effective Lorentz force with the gradient force in low 
energy solutions to the coupled Schroedinger equation. 
This cancellation is similar to that which occurs when 
appropriate external electric and magnetic fields “select” 
the velocity of a charged particle beam. This theoret¬ 
ical model illustrates how geometric magnetism could 
be harnessed to entertain similar capabilities for neutral 
systems. We then focus our discussion on neutral spin- 
1/2 systems subjected to an external magnetic held. We 
show, with the appropriate external held conhguration, 
how such a system can exhibit topological Aharonov- 
Bohm[27] (AB) behavior. Though it has been discussed 
previously [m US] [Hj within the framework of the BO 
approximation, here we offer fully time dependent, cou¬ 
pled channel, wave packet solutions without employing 
assumptions based on adiabaticity. This capability al¬ 
lows us extend the analysis to cases in which the BO ap¬ 
proximation is no longer valid. We investigate how this 
systems transitions at low energies, and one that allows 
an Abelian description in which the AB gauge vector po¬ 
tential emerges, to that in which non-Abelian features 
arise. Finally, we explore the relationship of this system 
to a phenomenon that occurs in molecular systems that 
possess conical intersections [4l [201 [ 2 ^ [29] . Atomic units 
will be used throughout, unless otherwise indicated. 


II. THEORY 
A. Velocity sieve 

Consider the model Hamiltonian for an atom, or n- 
state, system 


H = - 


2m 


VR^+Had{R)- 


(5) 


U{R)HboU^R), where Hbo{R) is an n-dimensional di¬ 
agonal matrix with eigenvalues ei{R), and U{R) is a uni¬ 
tary matrix. In this discussion we consider the case where 
n = 2, and R — {x,y) and limit Had{R) to be time in¬ 
dependent. We require U{R) to be single-valued [50] and 
express it 

U = exp(—icr3(//2) exp(—zcr2f^/2) exp(i(T3(//2) (6) 

where Ui are the Pauli matrices, and (/, H are single¬ 
valued functions of the planar coordinates (x^y). 

We choose, as in Ref. m, 

y) = + tanh(/3a:)) 

(t){x,y) = LBoy (7) 

where L,Bo,l3 are constants. 

In solving the Schroedinger equation ihd'ip jdt = H ij; 
it is useful to expand ip in the basis of the adiabatic 
eigenstates of Had- In that description we arrive at the 
set of coupled, Schroedinger-like, equations for the multi¬ 
channel amplitudes that are minimally coupled to a pure, 
non-Abelian, gauge potential A{R) = iW{R)^U{R), 
as well as the diagonal scalar potential matrix Hbo 
whose entries are labeled by €i{R) and correspond to 
the Born-Oppenheimer energies of Had- We now invoke 
the Born-Oppenheimer (BO) approximation and project 
these coupled equations unto, open, ground-state BO am¬ 
plitude F{R). We get. 



FIG. 1: (Color online) Classical trajectories (dotted 
lines) obtained from solutions to Eq. (12 1 for various 
initial velocities. Colored spheres correspond to, time 
lapse, expectation values obtained from solutions of the 
coupled channel quantum mechanical problem 
generated by Hamiltonian Eq. (§. The black spheres 
pass undeflected with their initial velocities. 
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(V - lAp^^F{R) + Vbo{R)F{R) = EF{R). 


R is the position operator for the atom and Had{R) is 
the adiabatic Hamiltonian that generates its internal dy¬ 
namics. Suppose it can be decomposed into the form 


VBo{R) = e2{R) + b{R) 

HR) = X- Ai2 ■ Ai2 (8) 
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where b{R) 0 [3T] is related, but not equivalent, to the 
adiabatic [32] or Born-Huang correction|33j. Typically 
the BO approximated is justified [55] if, for total (colli¬ 
sion) energy E, the inequalities ei(-R) >> E > € 2 {R), 
Ae = ei — 62 » b{R) are satisfied. Ap is a vector 
potential and is obtained by projecting the non-Abelian 
gauge potential so that Ap = Tr PAP P = \g){g\, 
where \g) is the ground eigenstate of Hpo with eigen¬ 
value € 2 {R), and Ai 2 is the off-diagonal component of 
A. With this parameterization we obtain, for the ground 
adiabatic state, an effective curvature 


H = V X HAp = 

k ^ttBoLP sech^{Px) cos{^tanh{Px)). (9) 

Though A describes a pure non-Abelian gauge and so has 
vanishing curvature, Ap may have a non-trivial curva¬ 
ture, which for a non-degenerate ground state discussed 
here has the Abelian form ff = V x Ap. 

The classical limit for Eq. Q corresponds to a situa¬ 
tion in which the motion of the atom is governed by[5] 

=vxH- VVbo (10) 

where R{t) is the atom position coordinate. 

In addition to the conventional scalar gradient force 
— 'Wbo{R), (sometimes called the Hellmann-Feynman 
force), the atom experiences an effective velocity de¬ 
pendent Lorentz force. We argued[T5l [15] that such 
forces could be exploited to construct effective “magnetic 
lenses” for neutral atoms or neutrons. Here we under¬ 
score that observation by demonstrating that both the 
induced scalar and vector forces could be used in con¬ 
junction to develop novel capabilities for the manipula¬ 
tion and control of neutral particle beams. In particular, 
we describe a neutral particle velocity sieve that exploits 
both the velocity dependent force arising from geometric 
magnetism as well as the Hellmann-Feynman force. 

It is well know that for charged particles a velocity se¬ 
lector can be realized by choosing a suitable geometry 
in which a uniform magnetic field, H, induces a Lorentz 
force that cancels the gradient force produced by an elec¬ 
tric field, E, so that for a particle of charge q and veloc¬ 
ity u, vqE[ — qE. Pursuing this analogy we construct 
a Hamiltonian of the form given in Eq. ([^ where the 
entries of the BO eigenenergies are given by, 

Hbo = 0-3 (e(H) -I- A) - b{R) 
e{R)=voyH{x) H{x) = \H\ (II) 


and where A is chosen to be a sufficiently large energy 
gap so that the BO projection approximation into the 
ground state is justified. In definition Eq. (II) we have 
included a counter-term b{R) defined in EqTQ. In the 
BO approximation it’s presence cancels the adiabatic cor¬ 
rection generated by the off-diagonal components of A. 
With this choice, Eq. ([To]) becomes 


= ^H{x{t)) + y{t)vod^H{x{t)) 


=-^H{x{t))+VoH{x{t)). (12) 

In the asymptotic region x —> — oo all forces vanish 
and if we take the initial condition y = 0 ,x = vq, 
x(t) = vot + xo,y{t) = 0 is a solution to Eqs. (HI- 
For other impact parameters and velocities numerical so¬ 
lution of Eq. ( [l^ predict trajectories in which the in¬ 
coming particle is scattered. This behavior is illustrated 
in Figure ([^ where dashed lines represent classical tra¬ 
jectories superimposed on the BO potential surface (ig¬ 
noring the energy gap A). The black line represents the 
solution in which an atom with initial velocity vq propa¬ 
gates at constant velocity unimpeded. The blue and red 
lines correspond to initial velocities slightly larger and 
smaller, respectively, than vq. Atoms with these veloc¬ 
ities are clearly scattered by the effective Lorentz and 
gradient forces. 

Below we demonstrate that this behavior is shared 
with quantum evolution generated by Hamiltonian <§• 
The atom is initially in its ground state in the asymp¬ 
totic region and we propagate its amplitude using the 
time dependent method given in Refs. [IS1[IS]. The pack¬ 
ets are described by 



whose components obey the coupled Schroedinger equa¬ 
tions 

r)n 

ih^^=- — V^g + V^^f-{V + b{R))g (13) 

where 

H = (A -I- e(i?)) cos(0(a;)) 

Vi 2 = exp(-i(/>(y))(A + e(R)) sin(H(a;)). 

(14) 

In the remote past the packets have initial starting posi¬ 
tions outside the interaction region with null impact pa¬ 
rameters. In Figure ([^ the black spheres represent prop¬ 
agation of packets with initial velocity Ug. The centers 
of the packets track closely the classical trajectory. The 
blue and red spheres have initial velocities slightly greater 
and smaller, respectively, than vq. The quantum simu¬ 
lation illustrated in that figure clearly demonstrates the 
proposed velocity selection effect for this system. Though 
Hamiltonian ([^ is a simple, time independent, two-state 
(or qubit) system, it’s experimental realization may pose 
challenges. Below we introduce another, more familiar, 
two-state system a neutral spin 1/2 system (eg. atom, 
neutron) subjected to a static external magnetic field. 
In the geometry discussed below we show how geometric 
phase induced, Aharonov-Bohm like, interferometry can 
be realized by it. 





B. Geometric phase, Aharonov-Bohm, 
interferometry 
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FIG. 2: (Color online) Time lapse illustration of a pair 
of coherent Gaussian wave packets initially in the 
ground state. At r = 0 the pair of coherent packets are 
shown in the upper right and left sides of the figure. At 
a later time, the packets coalesce thus creating the 
interference pattern. 


Consider the external magnetic field 
B = B{p) ^ + Bok 


(15) 


where 4 >, p are the polar and radial coordinates in a cylin¬ 
drical coordinate system, and Bq is a constant, and we 
ignore motion in the z directions as it can be factored 
from the planar motion. If B{p) = X/p then Eq. (15) 


describes the field generated by a wire with current along 
the z-axis axis superimposed with that of a homogeneous 
magnetic field Bq k. The Hamiltonian of a neutral spin- 
1 /2 atom is then given by 


H = - ^ ^ Pb a- ■ B 

2 m 


(16) 


(alternatively, for a neutron is replaced by its mag¬ 
netic moment). We can re-express the internal Hamilto¬ 
nian Had = pbct ■ B as 


_ , Bo -i exp{-i(j))B{p) \ _ 

I * exp{icj))B{p) -Bo ) - 


UHboU^ 

where 


(17) 


„ f y/B^ + B^p) 0 1 

"'’““'“’I 0 -/ByrPw" ’ 


and 


U = exp(—itT3())/2) exp(zCTin(p)/2) exp(jtT3(/)/2) (19) 


. _ 1 / cosH(/9 ) —1 ie *'^sinH(p) \ 
^ 2p \smn{p) 1 —cosH(p) J 



and so. 


Ap = Tr PAP = 


1 — cosn(p) 
2 p 



^JbI + b^{p) + b{p) 

sin^il(p)/2 ^l'{p)\ 


( 20 ) 


( 21 ) 


The effective curvature for the ground adiabatic state is. 


H = V X HAp 


k-^sm{n{p))n'{p). 

2 p 


( 22 ) 


The above analysis is relevant in studies of the motion 
of cold atoms, that have a magnetic dipole moment, in 
the vicinity of current carrying wire (or nanotube) [33] . 
Here we focus on a special case, that in which B{p) is a 
constant Bp. For this case, according to Eq. ( |2^ , H = 0, 
and because Vbo{p) is also constant, and ignoring b{p) 
(which is very small in the region p ^ 0 traversed by the 
packets), the atom does not experience either a velocity 
dependent Lorentz, or scalar force. However, the induced 
vector potential Ap is not trivial and is given by 


Ap 



(23) 


where 


$ = 7r(l — cosH) = 7r(l — 


Bo 


Bl + Bj 


(24) 


Packet propagation is again described as in Eq. (13) but 
now 


V = pb Bo 

Vi 2 = -ipB exp{-i(j))Bp. 

(25) 

In order to demonstrate the proposed thesis we prop¬ 
agate two identical, coherent wave packets as shown in 
Figure]^ The packets are displaced from the origin and 
are allowed to propagate, having been given initial veloc¬ 
ities that allow them to coalesce, at time Tc, and form the 
interference pattern shown in that figure. In our simula¬ 
tion we have first set Bp = 0 so that the packets prop¬ 
agate freely. In that case the horizontal line that passes 
through the center where the two packets meet, the wave 
function has the analytic form (see Appendix A) 


with (j), p, the azimuthal angle and radial distance, in a 
cylindrical coordinate system, respectively and tan(H) = 
B{p)/Bo. Thus, 

A{R) = iU\R)VU{R) = ^A^ + pAp 


Sa'^k'^ exp(— 


o 2 j 2 2 

2a k Y) \ 
4 a'^fc^+ 77 Q ' 


-|- TrT^g 


cos^{kr]) 


(26) 


where rj is the horizontal coordinate, fc is a scaled 
wavenumber, rjo are the initial displacements of the pack¬ 
ets from the origin and a is the initial width of each 
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packet. This function is plotted by the red line in Figure 
([^. In panel (a) of that figure we superimpose the values, 
shown by the blue circles, obtained in our numerical sim¬ 
ulation. We find excellent agreement with the analytic 
result Eq. (26) and validates the numerical procedure 


used in this study. In panel (b) we plot the correspond 
interference pattern when B^^Bp has the value so that 
n = tt/2. We notice a distinct shift in the calculated in¬ 
terference pattern from that given by Eq. (26). However, 
a fit with the replacement in Eq. (26) 


cos^ (kr]) —)■ cos^ (kij + $/4) 


(27) 


provides an excellent approximation to the calculated 


ft A ft 
A t A 

WVVW 
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FIG. 3: (Color online) (a) Cross section of the 
interference pattern shown in Figure ( (^. Red dashed 
line represents the analytic form Eq. (26). Blue dots 


represent results from numerical solution for the case 
— 0. (b) Blue dots represent interference cross section 
for the case 12 = 7r/2. (c) Same as above for 12 = tt, (d) 
Same as above for 12 = 37r/2. 


data given by the simulation. 

We now consider the propagation of the two initial 
packets, at 2 = 0, having their vertical velocities reversed 
so that they propagate into the upper half plane and 
coalesce at a point that is the reflection of the coalescence 
point shown in Figure ([^. The resulting interference 
pattern is illustrated in Figure (|^. At the point d the 
pattern is again well described by Eq. (26) but with the 
replacement 


cos^{kf]) ^ cos^{kr] — ^/A) (28) 

Because Hamiltonian ( |16| ) is not invariant under reflec¬ 
tion about the ry (horizontal) axis it is not surprising that 
the interference pattern at point d differs from that at 
point b. It might not be as obvious that the difference is 
topological in nature. According to the discussion above 
the locations on the r] axis where local minima occur is 
given by 


Vm — 


rriTT ^ $ 
^ ^ 4fc 


(29) 


where the ± sign identifies the points on the horizontal 
lines passing through the points b and d respectively and 


[llllll 
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FIG. 4: (Color online) Fringe shifts produced by wave 
packet propagation of ground state solutions for the 
Hamiltonian given in Eq. (16) 


m is an integer. Therefore there is a displacement 

^ (30) 

shown by the right hand panel of Figure (|^, in the loca¬ 
tion of the relative minima between the upper and lower 
fringe patterns. It depends on the quantity 4>, which 
according to AB theory is given by 


dr ■ Aab 


(31) 


where C is a contour that encircles the path adcba in Fig¬ 
ure (|^ and Aab is given by Eq. (23). The connection with 
AB theory, and the topological nature of the fringe shift, 
becomes evident when we shift the packet paths, so that 
a displaced closed circuit a'd'c'b'a' no longer includes the 
fictitious flux tube located at the origin. In that case our 
simulations show that the difference in the fringes at the 
corresponding locations of b'^d' disappear, in harmony 
with the predictions of the gauge theory analysis. 

In the discussion above we presented the full quantum 
mechanical simulation of coupled equations (13) in order 


to derive the fringe patterns discussed above. However, 
in order to gain additional insight it is useful to appeal 
to semiclassical analysis of this system. Below we show 
how the results discussed above can be re-derived using 
a semiclassical framework in the adiabatic representation 
ofEqs. (p|. 


C. Semi-classical description of packet propagation 
in an Abelian gauge potential 


Consider the wave packet defined by Eq. 


(A1) and whose center, at r = 0, is located at point a in 
Eigure (|^. We need to predict the packet that grows out 
of it and whose evolution is determined by the coupled 
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Schroedinger Eq.(13). Our calculation show that, under 
the adiabatic condition A/fc^ >> 1, A = /tb J + Bj,, 


and in the adiabatic gauge, a good approximation for it 
at the time its center arrives at b is 


U{a,b)t{ji{b) (32) 

where i/'i is the free particle packet and the unitary op¬ 
erator U (a, b) is given by 


17(a, 6) « exp(i / dr ■ Ap). 
Jci 


(33) 


Here Ci represent a path integral along segment a—b that 
starts at a and ends at b and Ap is the gauge potential 
given by expression (23). 


Similarly, a packet initially centered at c at r = 0, 
translates along path c — b and arrives at 6 at Tc. It can 
be expressed 


Uic,b)Mb) 


(34) 


where ^”2 is defined in Eq. (A1). The coherent sum of 


these amplitudes at r = is then given by the expres¬ 
sion 


V'(&) = U{a,b)tlji{b) + U{c,b)tp 2 {b) = 

U{c, b){U-\c, b)Uia, 6)V'i(&) + Mb)) = 

U{c,b){U{b,c)Uia,b)Mb) +Mb)) = 
U{c,b){U{a,c)il)i+Mb)) (35) 


where we made use of the unitary property of U and the 
relation U{a,c) = U{b,c)U{a,b). Therefore 

\M)\'^ = \i^2{b) + U{a,c)Mb)\^ (36) 


Evaluating 

U{a,c)=exp{i / dr ■ Ap) = exp{—i — ) (37) 

J abc ^ 


and inserting this into Eq. (35) we obtain expression Eq. 

m. 


We obtain an analogous relation for the case where 
the momenta, along the ^ direction, of the initial wave 
packets at a, c are reversed so that at time Tc the packets 
meet at point d in Figure 4. Following the steps outlined 
above we find 


\i;'{d)\^ = \Md) + U'{a,c)Md)\^ (38) 

where tp'i are the corresponding free-particle packets 
whose ^ momenta are reversed, and 


U'{a,c) = exp(i / dr ■ Ap) = Uw U{a,c) (39) 
J adc 

where Uw is a Wilson loop integral 

Uw = exp(i j) dr ■ Ap) = exp(i$) (40) 

and the closed, counterclockwise, circuit encloses the ori¬ 
gin. 


According to Eq. 
therefore. 


(A4) 


and 


\^'id)f = \Mb) + UwU{a,c)Mb)f- (41) 

So the interference pattern |■^/;'(d)p at the top panel in 
Figure (|^, differs from that at the lower panel by a, 
gauge invariant, phase determined by the Wilson loop 
Uw, in harmony with the results obtained by the, fully 
quantal, numerical simulation. 



FIG. 5: (Color online) Interference fringes for energy 
defect A/fc^ = 0.1 and 11 = 27r/3. Panels (a),(b) are 
fringes (blue points) at locations b,d (bottom,top) in 
Figure (|^ of text. Panels (c),(d) show fringes for 
(excited) state selected measurements at the latter 
locations. The green vertical line is a reference line to 
aid the eye in comparing fringes between the top an 
bottom panels. The dashed red lines correspond to 
fringes predicted by the Abelian semiclassical theory 



FIG. 6: (Color online) Circuit C which encloses 
fictitious flux tube (red disk), is translated to new 
circuit C" that does not enclose it. 


D. Semi-classical description of packet propagation 
in a, pure, non-Abelian gauge potential 

If the collision energy, is much larger than the 

energy defect 2A, between the Zeeman split spin states, 
non-adiabatic transitions between those states can occur. 
Thus if the initial localized wavepacket, say at point a 
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FIG. 7: (Color online) Interference fringes for energy 
defect A/fc^ = O.I and 11 = 27r/3 for loop diagram C in 
Figure Panels (a),(b) show fringes (blue points) at 
locations b,d (bottom,top) in Figure (|^ )respectively. 
Panels (c),(d) show fringes for (excited) state selected 
measurements at points b,d (bottom,top) in Figure Q 
respectively. The green vertical line is a reference line to 
aid the eye in comparing fringes between the top an 
bottom panels. 


in Figure Q, describes a particle in the ground Zeeman 
level, it will not necessarily stay in that level as the packet 
evolves in time. 


In the discussion above we considered the adiabatic 
limit in which the ratio of energy defect to collision en¬ 
ergy is large i.e. A/fc^ >> 1. In that limit spin flipping 
transitions between ground and excited Zeeman levels 
are suppressed. We showed that, in this limit, the sin¬ 
gle channel (or Abelian) Schroedinger equation ^ with 
gauge potential Eq. (23) and Vbo = 0 accurately pre¬ 
dicts wavepacket dynamics and topological AB features. 
As the collision energy is cranked up so that 2 Ajk"^ < 1 
we anticipate that the single channel (Abelian) descrip¬ 
tion breaks down and non-Abelian features arise. 


In Figure ([^, panels (a), (b) we plot the interference 
patterns, corresponding to the top and bottom regions 
shown in Figure ([^ for the collision energy correspond¬ 
ing to A/k'^ = O.I. Though the collision energy is suffi¬ 
cient to cause Zeeman level transitions, our results sug¬ 
gest that many of the Abelian features persist. First 
we note that there is a phase shift between the inter¬ 
ference patterns, for the top and bottom regions respec¬ 
tively, that is nearly, but not exactly, predicted by the 
Abelian AB theory (which are shown in red in that fig¬ 
ure). Panels (c),(d) show the fringes, for the top and 
bottom regions respectively for state-dependent proba¬ 
bilities, in this case for excitation into the upper Zeeman 
level. A distinct phase shift in the fringe patterns for ex¬ 
citation is also seen, though its structure is not predicted 
by the Abelian theory. 

In order to investigate whether these features are topo¬ 
logical we translate the loop C, shown in Figure (^, into 
the loop C' and repeat the calculations described above. 
Because loop C" does not enclose the fictitious flux tube 


(shown by the red disk) classical AB theory suggests that 
the difference in fringe shifts, evident in Figure ([^, is 
null. Indeed, this is the case. In panels (a),(b), of Figure 
the probability interference are shown. The patterns 
for top and bottom regions (points d,b respectively) line 
up and no fringe shifts are evident. Panels (b),(d) of 
that figure show the corresponding interference patterns 
for the excited Zeeman level probabilities. Though fringe 
differences are negligible, we note a strong suppression of 
the latter (when compared to that for loop C shown in 
panels (b),(d) in Figure ([^ ). The suppression of excita¬ 
tion for loop C is clearly a non-Abelian (or multichannel) 
feature. 

Finally, we consider the extreme non-adiabatic regime. 
In it the ratio Ajk^ —)■ 0 and we can again employ semi- 
classical methods [33] to predict the fringe patterns that 
are generated by our fully quantal simulations and which 
are shown below. In Figure (|^ we repeat the calculations 
for propagation along loop C in Figure 0 for the values 
A/k"^ = 0.007 and = 2tt/2>. Unlike the case for the 
adiabatic and near adiabatic regimes, in which topologi¬ 
cal fringe shifts arise, panels (a), (b) of Figure ^ clearly 
demonstrate absence of the topological fringe shift. This 
behavior can be explained using semiclassical methods. 
For, in that description [34] the total probability ampli¬ 
tude at point (b) in that grows out of wave packet 



FIG. 8: (Color online) Interference fringes for energy 
defect A/fc^ = 0.007 and U = 27r/3. Panels (a,b) fringes 
(blue points) at locations b,d (bottom,top) in Figure (|^ 
) respectively. Panels (c,d) fringes for (ground) state 
selected measurements at points b,d (bottom,top) in 
Figure (|^) respectively. The green vertical line is a 
reference line to aid the eye in comparing fringes 
between the top an bottom panels. The dashed red 
lines correspond to fringes predicted by the 
semiclassical Abelian theory 


ipiia) is approximated by the expression 


Pexp(z / dr ■ A)'tjji{b) 


(42) 


where P is a path-ordered integral, or Wilson line, along 
segment a — b, and A is the non-Abelian, pure, gauge 
potential (20). Repeating the argument outlined above 






























for the fringe shift in the adiabatic limit, we now find 
that the shift depends on the Wilson loop integral 


Uw = P exp{—i dr ■ A) 


(43) 


where A is given by Eq.(20). Because A describes a 


pure gauge we find that Uw = 1 (see Appendix C) and 
so, unlike the case in the adiabatic regime, a fringe shift 
between the interference patterns at 6, d in Q does not 
manifest. Interestingly, this is no longer true if we per¬ 
form state dependent measurements at locations b, d in 
that figure. In panels (c), (d) of Figure (|^ we plot the cal¬ 
culated interference patterns at those points for a state 
selective measurement (in this case, the ground state). 
Those fringe shifts are, again, accurately predicted by 
the Abelian theory as illustrated by the red lines in those 
panels. Furthermore, this shift is also topological, in 
that the shifts vanish for loops that do not enclose the 
fictitious flux tube. A detailed discussion of the origin 
and implications of this observation will be presented 
elsewhere [34] . 



FIG. 9: (Color online) Standard, single loop, AB 
inteferometry setup. Two coherent packets at origin a 
split and propagate to the mirrors at points b, d. The 
mirrors deflect the packets so they recombine at point c 
where the measurements are made. 


above we find measurements at screen c are now propor¬ 
tional to the loop integral 


P exp(z 


dr • A) = 1. 


(46) 


E. AB inteferometry for a single loop 

Our discussions addressed AB inteferometry for se¬ 
tups in which interference patterns are compared follow¬ 
ing two independent open-loop measurements |55] . In the 
classical single loop AB setup, an interference pattern is 
observed at a single screen (point c) as shown in Figure 
In it, a wavepacket coherently splits at the origin, 
point a, and propagates toward the mirrors at points b,d 
respectively, the packets are deflected and allowed to re¬ 
combine at point c where a measurement is taken. In the 
calculations described above, and in the adiabatic limit 
A/fc^ >> 1, we found that (i) the wave packets propa¬ 
gate as a free particle, (ii) in the journey along curve C 
the packets acquire, in addition to the standard dynam¬ 
ical the phase factor, the phase 


exp(i / dr ■ Ap). 

Jc 


(44) 


Using Eq. (44) for the paths a — b — c and a — d — c we 


find that measurements at the screen will be a function 
of the path integral, along the closed loop, 


exp(j dr ■ Ap). 


(45) 


If that loop encloses the fictitious flux tube then its value 
is exp(i<I>), otherwise it has unit value. Therefore this 
setup exhibits topological properties consistent with AB 
theory. In the extreme non-adiabatic limit A/fc^ << 1 
our calculations again demonstrate the validity of prop¬ 
erties (z), (zz), with the exception that “free particle” evo¬ 
lution is that of a two-component wavefunction and the 
phase factor multiplying it is a multi-channel unitary ma¬ 
trix Eq. (42). Therefore, repeating the analysis given 


That is, regardless of the loop geometry, topological 
fringe patterns do not arise and there is no topological AB 
shift. In summary, our calculations demonstrate that, in 
the adiabatic limit, solutions generated by the time evo¬ 
lution operator Eq. ([l^ with the external field configu¬ 


ration Eq.(15l reproduces the standard AB fringe shift. 
In the limit A/fc —>■ 0, non-adiabatic transitions conspire 
to wash out topological AB fringes. This conclusion is 
consistent with that given in Ref. [TSj. Neverthless, we 
find here that if spin-state dependent measurements are 
made during traversal of the circuit, topological shifts 
persist. This counter-intuitive observation will be dis¬ 
cussed in more detail elsewhere 1341. 


III. AB SCATTERING 

In the previous section we presented an outline for elic¬ 
iting AB-like topological effects, for neutral spin 1/2 par¬ 
ticles, in an interferometry setup. In standard treatments 
e.g. [36ll37], the AB effect is discussed in the context of 
a scattering scenario and so it is instructive to analyze 
the dynamics generated by Hamiltonian Eq.(16) in this 


framework. In Ref. |18] a Born-Oppenheimer approxima¬ 
tion was used to obtain an equation for the ground state 
amplitude (in a field theoretic analog of this system), and 
in which, the vector potential Aab emerges. Standard 
time independent methods were applied to demonstrate 
AB -like scattering. A BO approximation was also used 
in Ref.[T2] in order to suggest that AB -like scattering 
arises in systems described by Hamiltonian (16). 


Here we apply time dependent wave propagation in 
order to demonstrate AB scattering. However, there are 
two obstacles that prevent us from simulating pure AB 
scattering in the Frauenhofer region. The latter requires 
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propagation of a wave packet from the distant remote 
past, to the distant future. This is numerically unten¬ 
able. In addition, the wave packet traverses a region near 
the origin where singular couplings between the ground 
an excited adiabatic states arise[31]. We by-pass these 
difficulties by placing an impenetrable hard-cylinder sur¬ 
rounding the origin and introducing a diagonal counter 
term hi{R) = 


Aji, as shown in Eq. (13). 


At initial time tg we introduce a finite slab-like wave 
packet shown by the illustration in Figure It pro¬ 
ceeds with a mean initial velocity toward the impenetra¬ 
ble cylinder represented in that illustration by the orange 
colored disk. The wave packet is propagated numerically 
solving coupled equations (13). The packet mean mo¬ 
mentum is chosen so that the inequality A/k^ >> 1 is 
satisfied. At a later time ti Figure illustrates how 
the packet diffracts around the cylinder. Subsequently, 
at t 2 , the forward scattering probability density reaches 
a detector shown by the red dashed line in that figure. 

In Figure [TT] we present a high resolution plot for the 
imaginary part of the wave packet amplitude, at time 
t2, for the two cases where fi = 1/2 (left panel) and 
n = 0 (right panel) respectively. With = 0 Eqs. 
(13) decouple,and the ground state amplitude satisfies 


the free particle Schroedinger equation in the presence of 
a impenetrable cylinder centered at the origin. The right 
hand panel of Figure [TT] illustrates both the reflected and 
transmitted components, including interference between 
the two, of the packet as it is scattered by the cylinder. 
In the left hand panel we plot the same amplitude except 


we set the parameter = 1/2 in Eq. (14). It corresponds 
to the case where $ = tt in Eq. (23). The plot is almost 


identical to the one obtained for 11 = 0, except for the 
striking phase dislocation[3S], running along a line bi¬ 
secting the cylinder in the forward scattering direction. 
This phenomenon is well know in pure AB scattering|36j 
and, for the case $ = tt, is a manifestation of a nodal line 
for the amplitude along the phase dislocation line. Berry 
et al.[35] argued that phase dislocation in pure AB scat¬ 
tering, though not a physical observable, is a topological 
invariant. In Figure (12) we plot the probability ampli¬ 


tude at the observation panel at time ^2, shown by the 
line in Fig. (10), for both cases H = 0, the blue dashed 


line, and 11 = 1/2 by the solid red line. Both probabil¬ 
ity distributions are similar except along the line near 
ry = 0 in that figure. Both cases exhibit the “shadow” 
cast in the forward direction by the cylinder, the fl = 0 
plot show a small enhancement directly ” downwind”, at 
77 = 0, from the cylinder. It corresponds to the ID analog 
of the celebrated Poisson spot [35113S] that occurs when a 
wavefront diffracts about a circular obstacle. The red line 


in Fig. (12), which corresponds to the case where $ = tt 


in AB scattering and which shows a strong suppression 
of this spot. 

Finally, we compare the results shown in Fig. 


( 12 ) 


with those obtained in a time independent description. 
In appendix B, we construct the wave functions used in 
a standard scattering theory scenario. It obeys boundary 


conditions that correspond to an incoming packet in the 
asymptotic region in the remote past that subsequently 
scatters off the cylinder. Typically this wave function is 
used to derive the scattering amplitude in the Frauen- 
hofer region, but here we use it to find the, time inde¬ 
pendent, probability amplitude at the observation screen 
shown in Fig. (10). In Figure (12) the dashed blue line 


represents the case 4) = 0 (i.e. no AB flux tube), and 
the red line by the $ = tt case. The latter corresponds 
to maximal AB scattering. As in Fig (12), differences in 


the probability distributions are most evident near the 
ry « 0 line. Because the incident amplitude extends from 
—00 < 7y < 00, Fresnel type interference patters at larger 
Tj manifest here and which are not present in Fig. (12). 


However, there is qualitative agreement between the two 
descriptions near the forward direction, in that a strong 
suppression of the Poisson spot is evident when the AB 
flux tube has the value 4) = tt. 

In summary we have demonstrated that solutions of 
Eq. (13) do exhibit the features that arise in standard 
AB scattering for the collision energy range A/k^ >> 1. 



FIG. 10: (Color online) Time series plot of the 
probability density for an initial wavepacket at tg that 
propagates toward the impenetrable cylinder (orange 
disk) and is scatter'd by it. The red dashed line 
represents a detection screen. 


IV. RELATION TO THE MOLECULAR AB 
EFFECT 


A gauge potential equivalent to Eq. (23) (for $ = 
tt) was introduced by Mead and Truhlar JT order to 
describe molecular dynamics near a conical intersection. 
In this section we explore the relationship between that 
phenomenon, often called the molecular Aharonov-Bohm 
(or MAB) effect, and AB-like scattering described in the 
previous section. It has long been taken for granted, in 
the molecular physics community, that degeneracy in the 
form of a conical intersections is an essential requireiuent 
for topological effects induced by a vector potential to 
arise. However, no degeneracy in the physical parameter 
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FIG. 11: (Color online) Plots of the imaginary part of 
the ground state amplitude at time t 2 shown in Figure 
(10). Right panel corresponds to case $ = 0, scattering 
by an impenetrable cylinder. The left panel corresponds 
to the case d) = tt, maximal AB scattering by an 
impenetrable cylinder. 





FIG. 12: (Color online) Dashed blue line is a plot of 
the time dependent probability density, at ^2 and for 
$ = 0, along the 77 axis corresponding to the red dashed 
line shown in Figure (10). The red solid line 
corresponds to the case $ = tt, maximal AB scattering 
by an impenetrable cylinder. 


space R is evident in the BO Hamiltonian defined in Eq. 
(18) as Bq, Bp are taken to be constant. This dichotomy 
presents a certain amount of cognitive dissonance and so 
deserves closer examination. 


Consider a tri-atomic system that possesses a conical 
intersection at the origin of a parameter space that is 
spanned by a set of nuclear internal coordinates x,y. 
Typically they represent various linear combinations of 
the squares of internuclear distances between the three 
nuclei]!] in a planar configuration. In this coordinate 
system the azimuthal angle (j) is called the pseudorota¬ 
tion and p = + ip' measures distortions from an 

equilateral triangle configuration of nuclei. We describe 
the system by an amplitude 'ip{x,y,r) where r are elec¬ 
tronic coordinates. If the electronic, or fast, coordinates 
are integrated out one can approximate the adiabatic, or 
electronic, Hamiltonian as a truncated two-dimensional 





FIG. 13: (Color online) Dashed blue line is a plot of 
along the y axis, corresponding to the dete ction 
screen shown in figure (10), for solution to Eq. (Bl) 
with $ = 0. It corresponds to time independent 
scattering from an impenetrable cylinder at the origin. 

The red solid line corresponds to the case $ = tt, 
maximal AB scattering by an impenetrable cylinder. 


Hilbert space operator, which in the vicinity of the inter¬ 
section is given by|55] 




The eigenvalues of Had are + y'^ and correspond 

to first excited and ground states, respectively, of the 
electronic Hamiltonian. We approximate the vibronic ki¬ 
netic energy operator Hke = + where p, 

is a reduced atomic mass. Therefore, 


2p V dx"^ 



+ Had 


(48) 


and has the form given by Eq. ([^. 

Because Had is real, Longuet-Higgins and Herzberg[29] 
constrained its eigenstates to be real-valued and found 




where |$g) is the ground adiabatic electronic state. They 
noted that it is multivalued, as its value changes sign in 
traversing a circuit from (p = 0 to (j) = 2%. The total sys¬ 
tem amplitude tp must be single valued and so in a Born- 
Oppenheimer approximation in which tp = F(x,y)\^g), 
the vibronic amplitude F{x,y) must undergo a compen¬ 
sating sign change. That argument was used by Mead 
and Truhlar to invoke the minimal coupling of the vi¬ 
bronic motion, in Eq. (48), with a vector potential given 
by 


AmAB 


2p 


(50) 
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In order to relate this result with the analysis in the pre¬ 
vious section we offer a different tack. Instead of con¬ 
straining the eigenstates of Eq. ( [47|) t o be real we do not 
impose phase restrictions on them|6j, but we do require 
them to be single valued in parameter space x, y. We 
find that 

Had = UaWHBoU^cW 

_ ( + y‘^ 0 


Hbo = 


0 — + y"^ 


(51) 


where 


Uc{(t>) = exp(-i(T2(/)/2) exp(-icr3((i/2) = 



(52) 


Unlike the operator U{(p) given in Eq. (49), which under¬ 
goes a sign change as (j) ranges from 0 to 27r, Uc{(j)) is sin¬ 
gle valued for all </>, excluding the origin, i.e. 17((/)-|-27r) = 
U{(j))- Repeating the derivations outlined in the previous 
sections we arrive at the set of coupled equations for the 
vibronic amplitudes in the adiabatic picture. 


2m V 


F{R) + Vbo{R)F{R) = EF{R). 

(53) 


E(R) is a column vector whose entries are the ground 
and excited state vibronic amplitudes, Vbo{R) is the BO 
eigen-energy matrix and A is the non-Abelian gauge po¬ 
tential 


1 


A = iUlVUc= — > . 


-1 


-iexp(-i<()) j 


It is identical to expression ( [20| ) for the case 0 = Sir/2. 
This equivalence is a consequence of the fact that [/(</>), 
and Uc are related, as shown in Appendix C, by a con¬ 
stant unitary matrix. A Born-Oppenheimer projection 
of Eq. (1^ leads to Eq. ([sl) where Ap becomes the 


Mead-Truhlar gauge potential Amab- 

U{R) defined in Eq. (19) diagonalizes the adiabatic 
Hamiltonian pbct ■ B, which for the B configuration Eq. 
(15) does not possess degeneracy in R space. Uc diago¬ 


nalizes Hamiltonian (47) which exhibits a conical inter¬ 
section at the origin of R space. Thus the emergence 
of the non-trivial gauge structures follow from the prop¬ 
erties of the unitary operators that diagonalize the re¬ 
spective adiabatic Hamiltonians, and is insensitive to the 
fact that Hbo may possess degeneracy. In Appendix C 
we show how both U (i?) and Uc can be constructed, via 
a Wilson line integral, from their respective gauge fields. 

Einally, we explore the relationship of this analysis 
with that given in discussions of Berry’s phase [1]. The 
adiabatic Hamiltonian Had = given in Eq. ( |16[ ), 

shares the structure of Berry’s model. In the latter B 


t / 
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EIG. 14: (Color online). The left panel illustrates the 
vortex configuration, projected on the xy plane, of the 
external magnetic field given by Eq. (15), and the red 


dashed lines are paths in R space. $ represents the 
effective flux tube pointing into the page. The right 
hand panel shows field lines emanating from the 
effective Dirac monopole located at the origin. The 
dashed red lines represent loops in the parameter space 
of the monopole, and are maps of the corresponding 
loops in R space. The blue background plane bisects 
the monopole. 


is taken as a classical parameter that can be arbitrarily 
varied. Under adiabatic conditions, i.e. slow variation 
of B, a spin 1/2 system accrues a phase that is deter¬ 
mined by the vector potential, or connection, of a mag¬ 
netic monopole (see the right panel of Figure 14 ). The 


gauge field of the monopole differs from that given in Eq. 
(23) which describes a magnetic flux tube extended along 


the z axis. In our discussion R is a fixed external mag¬ 
netic field, shown by the vortex lines in the left panel of 
Figure and gauge field Eqs. ( 20|23 ) follows from it by 
considering variation of parameter R instead of B. The 
gauge structure induced by the former, realized in molec¬ 
ular systems with conical intersections, does not lead to 
a Lorentz force. However, effective magnetic monopoles 
in R space are realized in diatoms [5l l6l [2T] . They lead to 
effective Lorentz forces, (i.e. geometric magnetism), in 
those systems[ni [H]. In addition, i? is a quantum vari¬ 
able and so true adiabaticity is ill defined. Nevertheless, 
if we demote i? to a c-number we can consider loops in R 
space and relate them to corresponding loops in B space. 
This is illustrated in Figure (14). In the left panel of this 


figure we show a loop Ci which encloses the flux tube, of 
magnitue $, at the origin. In B space it corresponds to 
a loop that subtends a solid angle U, from the location 
of the monopole, so that <I> = 7r(l — cosfl)/2. A loop 
C 2 that does not enclose the flux tube corresponds to a 
loop in B space that lies in a plane (shown by the blue 
background in the right panel of that figure) that bisects 
the monopole. 
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V. SUMMARY AND DISCUSSION 


particle wave packet (normalized to the value 2) 


These examples illustrate how induced geometric phase 
gauge potentials, that lead to holonomies and effective 
“magnetic” forces, arise in simple quantum systems with¬ 
out resorting to assumptions that are based on BO - type 
approximations. We introduced a two-level Hamiltonian 
which mimics the behavior of a charged particle that is 
subjected to both electric and magnetic fields in such a 
way that enable a velocity selector. In another exam¬ 
ple, that of a neutral spin 1/2 system subjected to an 
external magnetic field, we showed how AB-like effects 
arise in low energy solutions to the coupled Schroedinger 
equations. We introduced and discussed both an inter¬ 
ferometry setup, as well as a scattering scenario, in which 
such effects arise. In addition we showed how the system 
transitions from that described by an Abelian gauge po¬ 
tential to that in which non-Abelian effects manifest at 
higher collision energies. 

Dirac |2] noted that if a solution ijj to the Schroedinger 
equation is multiplied by a phase factor so that —>■ 

ip'iR) = exp{ia{R)) '/’(H), and since we can always rep¬ 
resent a{R) = dr ■ A(r) this substitution leads to the 
replacement V')/ —>■ — iAijj'. However, conventional 

quantum mechanics demands that ijj'{R) be single-valued 
at all R, and this condition constrains the gauge poten¬ 
tial A to a trivial, pure, gauge which does not lead to 
new physics. Dirac suggested that minimal coupling of 
a charged particle with a non-trivial gauge field follows 
from a non-integrable factor a(i2, C), a functional that 
depends on the path C. There have been explorations 
in that direction[40] but difficulties in enforcing single¬ 
valuedness has limited the utility of this approach. 

The examples provided above illustrate how an as¬ 
pect of Dirac’s program is realized. Let '0 be a multi- 
component amplitude and consider the non-Abelian ver¬ 
sion of Dirac’s substitution, i.e. 0 —>• 0' = [/(i2)0, where 
U{R) is a single-valued, differentiable unitary matrix op¬ 
erator. In m and in Appendix C we showed how U(R) 
can be expressed as a Wilson line of a non-Abelian, pure, 
gauge potential A. Because A is a pure gauge, issues 
involving multi-valuedness in U{R) and hence in 0' do 
not arise. Nevertheless, as our examples demonstrate, if 
the gauge symmetry is broken by energy gaps (i.e. non¬ 
degeneracy of HBO ) non-trivial gauge fields that lead to 
gauge forces and/or topological holonomies can emerge 
in low energy solutions to the Schroedinger equation. 


Appendix A: Packet dynamics 


At lower collision energies, in which the excited Zee- 
man level is closed, we consider the propagation of a co¬ 
herent wave packet that is initially localized, as shown in 
Figure (i). Using the dimensionless coordinates defined 
in [TS], ^ = x/L,rj = y / L,T = 2m 


= 0l(^,r7,T) -h02(C,?7,t) 

0i = y y dkidk2^{ki — kxi)^{k2 — ky. ) X 


exp(—ir(fcj -I- fc|)) exp(—irA) x 
exp{iki{^ -^i)) exp{ik 2 {r] -Vi))- 

(Al) 


With the choice 


^{k) = -^ exp{—a^k^) 

(271^) 

Vi x) is a solution to the equation 



5^0, 

de 




+ A0i 


(A2) 


(A3) 


and describes, at r = 0, a Gaussian wavepacket of 
width a localized at i^i,Vi) with momenta kxi^ky-. A 
is the energy defect A in units of By linearity, 

the coherent sum 0(^, ? 7 ,t) is also a possible packet. 
Here we chose Ci = '?2 = 0,Vi = —V 01 V 2 = Voj 
kxi = kx 2 = k, ky-^ = —ky^ = k. At time Tc = ??o/2fc the 
two packets coalesce and form the interference pattern 
illustrated in Figure The center of the merged packets 
at Tc passes through the line ^ = Vo- On it 


0(f?O,?y,Tc) 


dkidk 2 $(fci — k)^{k 2 — k) x 


exp(—tTc(fci -I- ^2)) exp(iA:i?7o) exp(—zrA) x 
exp{ik2{v + Vo)) + exp{-ik2iv - Vo)) - 


(A4) 


Integration of Eq. (A4) yields 

I0(?7o, ??,Tc)P = 


4 a ‘* fc ^ 7 r -|- ttvq 


cos^{kv)- 


(A5) 


Appendix B: AB scattering wave fnnction in the 
presence of an impenetrable cylinder 

Scattering by an impenetrable cylinder, in the pres¬ 
ence of a AB flux tube, has been extensively reviewed in 
the literature [MIMlISZl- Below we provide an analysis 
following along lines given in Ref. [36]. We consider the 
Schroedinger equation, 

(V - - U(r)i, = 0 

A = (Bl) 

where A is a vector potential, and U{r) is a short range 
potential, wich we take to be that of an impenetrable 
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cylinder whose axis is perpendicular to the incoming flux. 
Expressing the wave function 

OO 

tp{r,(j))= ^ Cm exp{im(j)) Rm{r) (B2) 

m— — oo 

we obtain the radial equation, 


^ , R'ir) (to - a)2 

R (r) H- - - ^ - R{r) 

(k^ - U(r))R(r) = 0 


(B3) 


where a = ^/2 'k. Its solutions are linear combinations 
of 

'^Im—a| -^Im —al (^^) ■ 


For the repulsive core we require R{a) = 0 and therefore 
the radial solution has the form 

( 1 ) 


R{r) = ciJ\m-a\{kr) + C 2 Wj_Akr) (B4) 


where 


C2 '^|m—a|(^^) 

Cl .(ka)' 

|m— q;|'' ■' 


(B5) 


We define the AB amplitude ([25) 

m—oo 

i)ab{r,4>)= ^ (-z)l”""“lj|„_„|(fcr)exp(iTO(/)), 

m=—OO 

(B6) 

and so choose ci = (—Therefore, according to 
Eq. (1B5| 


C2 = -(-*) 




I [ka) 
\m—<y.\ ^ > 

and 

(j)) = ipabir, (j)) - -tpRir, (j)) 


(B7) 


\m—a\ V / 


(B8) 


For the value a = 1/2, the so called maximal AB scatter¬ 
ing case, ipabir, </), is given by the closed-form expression 


'!/)af,(r, 0) = exp(i(//2) exp(—i/cz) exp(—iTr/d) x 


where the ± refers to the regions 0 < /) < tt and — tt < 
4 > <0 respectively. 

Appendix C: Representation of a unitary matrix by 
a pure gauge field 

Consider a unitary matrix U(R) that is continuous, 
single valued, and has a derivative defined at each point 
R, excluding the origin, so that 


A{R) =iU\R)VU{R) 


(Cl) 


is a pure gauge potential[3]. Suppose it has the property 


Pexp(f ® dR • A) = 1 


(C2) 


for any closed path in R space. Can we then represent 
U (R) as a Wilson line of the corresponding gauge poten¬ 
tial ? Here we show that this is the case for the operators 
U{R) given in Eq.(19), and Uc defined in Eq. (52). We 
define a path C so that R{t) = Rq cos{ujt)i + Rq sm{ojt)j 
on it, and traverses a loop of radius Rq in the counter¬ 
clockwise direction. Consider the (Wilson) line integral 
along this path from (pQ to (p 


W{(p,(po,C) = Pexp{i ( dR-A[R)) = 

Jc 

W{t,to) =Texp{i f A{t)) 

Jtn 


(C3) 


where, (p/uj = t,(po/u! = to, T is the Dyson time -ordering 
operator. 


A{t)^^-A{Rit)) = 

—a;sin^(^) sin(r2) 


— Ize'^wsin)!!) 


wsin^d) 


(C4) 


and we used Eq. (20). Noting that 


dW{t) 


dt 


= iA{t)W{t) 


(C5) 


exp(=Fz7r/4)Ar/(exp(=Fz37r/4)-\/2/cr cos((//2)) (B9) we integrate to get 


T/f/ZX, X, i TTU^'\TTirh\ ( cos^ (§)-I-sin^ (§) e -e * 2 “ sin(^-^) sin(D) \ 

irw,0„) = ci'ww„)=l co.qi)+»nq|)e-«.-« )■ 


(C6) 
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W(4>,4>o) obeys the group property 

for cj) > cj)i > (fiQ, and since W is unitary we have 

exp(i f A{t)) 

Jto 

A{t) = ^-A{R{t)) (C8) 

where T is the anti-time ordering operator, and R{t) de¬ 
fines a circular path running along the clockwise direc¬ 
tion. According to Eq. ( |C6[ ) 

= (C9) 


= Texp^—i f A(t)) = T 

Jto 


shown that U{R) can be defined in terms of a line inte¬ 
gral along the path shown by the red dashed arc in Fig. 
(15). Point a is along the line (pQ (not including the ori¬ 
gin) and R is some arbitrary point (not on the origin). 
However we can deform the path as shown by the solid 
l ine segments in that figure. By the group property Eq. 
(|C7|) we can equate 


W{R,cl>o) = W{R,d)W{d,c)W{c,b)W{b,a) (CIO) 


where we used W{c, b) = W{R, d) = 1, and demonstrate 
that U{R) can be expressed in terms of the Wilson line, 
for any right handed path, provided that it does not pass 
through the origin. We stress, again, such a representa¬ 
tion is not possible for the Abelian projection 

Ap = Tr PAP = (cil) 

2p 

as its Wilson loop integral, for any loop surrounding the 
origin, is not unity. As a consequence 


where U{(f>) is given in Eq. (19). Therefore we have 



FIG. 15: (Color online). Integration paths for Wilson 
line integral, where point a lies on a ray p = pQ and R 
lies on the ray 4>. 


Up = exp{i / dR-Ap) 
Jc 


is multivalued. In the same way we also show that oper¬ 
ator Ucicp), defined in Eq. ( |^ is given by 

(C12) 


UM = W\<P,0) 


where we have set H = 37r/2 in the expression for W. 
Comparing Eq. (C12), setting H = 37r/2, with Eq. (C9) 
we find that 


U{P) = UiO)Uo{cP) 

U{0) = exp(37r*/2 CTi). (C13) 
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